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Abstract 

We study the loop induced rare decays Bs{B^) "fuv within the standard model. 
Both constituent quark model and pole model are used and the branching ratios turn 
out to be of the orders of 10^*^ for ^vv and 10~^ for B^ 71/P. These processes 

can be served to determine the decay constants of Bg and B^- 
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1 Introduction 



Rare decays of bottom quark ^ can be served as tests of the standard model (SM). As the 
observed process b ^ s'-f which induced B K*^ [|| and B Xs^ [H at the hadronic 
level, b are also studied extensively in the past both at inclusive and at exclusive 

levels. It was recognized that the SM prediction for the branching ratio of i? ^ Xsiyiy, which 
is ~ 5 X 10~^|^, is within one order beyond the present detectivity and this process might be 
observed soon0. Furthermore, the corresponding exclusive decay B — » K*vi' is important 
as an input to determine V^h to an accuracy of 10% 0. 

Pure leptonic decays of heavy pseudoscalar mesons into light lepton pairs are helicity 
suppressed, their branching ratios are 

B{Bs yuV") = 1.8 X 10~^ 

B{B, e+e-) = 4.2 x 10-^^ (1) 

which make it difficult to determine fs^ from these processes. For Bd the situation is even 
worse due to the small CKM mixing angles. Although the processes Bs{Bd) t^t^ do 
not suffer from this suppression mechanism and the branching ratio is about 10^^ in SM, 
it is hard to be detected at future i?-factory where the efficiency is not better than 10~^. 
The processes Bs{Bd) — > 7^^!^, whose branching ratios depend quadrically on fBsifs), can 
be taken as a possible alternate to determine the decay constants. For detections of these 
processes the method of searching for the "missing mass" can be used, which is similar to 
the case of i? — > X^z/z/ 

In the present work, we study the processes Bs{Bd) — > "^vv within the SM. We will 
analyze Bs{Bd) —>■ 7z/z/ in the following: In section 2 the relevant effective Hamiltonian will 
be given in SM. Constituent quark model and pole model will be used in section 3 to give 
the predictions. Finally, section 4 contains some brief discussion. 
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2 Effective Hamiltonian 



Let us start with the quark level process b qyv, with q = s or d. The Feynman diagrams 
are displayed in Fig.l. Both box and Z penguin diagrams contribute to this process. The 
resulting effective Hamiltonian is left handed in SM ||^, |^ : 

n = C{q-f,PLb){uYPLu), (2) 

with Pl = {1 — 75)/2. The coefficient C is 

where x = mf/m^. For simplicity, we neglect the QCD correction to this coefficient, whose 
effects are within 2% if appropriate renormalization point is chosen 0. 

Since the neutrinos are massless in minimal SM, the processes Bq{hq) ^ z/z/ Q are forbid- 
den by helicity. If an additional photon line is attached to any of the charged lines in Fig.l, 
the situation will be different: no helicity suppression exists any more. However, when the 
photon line is attached to one of the internal charged lines, there will be a suppression factor 
of ml/m^r in the Wilson coefficient compared with the ones for h qvv. The reason is that 
these effective operators are now dimension-8 instead of dimension-6. On the other hand, 
the diagrams in Fig. 1 with photon line connected to one of the external (bottom or strange 
quark) lines, whose effective Hamiltonian turn out to be 



n = -lcq 

o 



will be the dominant contributions to the decay Bq{hq) — » -^uu. 



h (uj^Pli^), (4) 



3 Model calculations 

The effective Hamiltonian given in is not enough to analyze the processes Bq ^vv. In 
addition, models are needed to do the calculations at the hadronic level. First we use a simple 
"'^We denote mesons with the quark content {bq) as Bq for convenience. 
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constituent quark model (see, for example [0). In this model both of the (anti)quarks are 
treated non-relativistically with the same velocity of the hadron. Thus (anti) quarks masses 
are now constituent quark masses of the order of several hundred MeV. We use further the 
interpolating field technique [0] which relates all the hadronic matrix elements of the present 
case to the decay constants of the mesons. The amplitude for Bq ^vv decay is: 



^ eCfB.mB. 
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l2[pBq-P-t) nib J 'II nib J 



(5) 

The first term inside the bracket comes from the CP-odd part of Bq, while the second term 
results from the CP-even part of Bq. Since mb » niq {q = d,s), we can safely neglect the 
term of l/m^. Then after squaring the amplitude, both CP-odd and CP-even parts give 
the same contribution to the decay rate. Performing the phase space integration over one 
of the two Dalitz variables, and summing over the three generation of neutrinos, we get the 
differential decay width versus the invariant mass: 

^Ba - ^vu)^vu ■ (6) 



(im^p 3(487r)2mB^m2 
The decay width is: 

(1447r)2m2 ' ^ ' 

Using a = 1/132, = 0.51 GeV, nit = 176 GeV and \VtbVt*\ = 0.04, we get 

r(S, ^ ^uu) = 9.5 X 10-21 X (^^) GeV. (8) 



If the lifetime is taken as t{Bs) = 1.34 x 10 |I0|, the branching ratio is found to be 
1.9 X 10^^. The differential decay rate versus m^p is displayed in Fig. 2 (solid line). 

For Bd meson decay, we take = 0.01 and nid = 0.35 GeV. The decay width is 

then 



T{B, ^ ^ujy) = 1.2 X 10-'' x ( 1 GeV, (9) 
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which corresponds to the branching ratio 2.6 x 10 ^, if uses of T{Bd) = 1.50 x 10 ^^silTOl and 
fs. = 0.2 GeV are made. 



We also apply the pole model to calculate the Bq meson decay. The main contribution 
here is the radiative transition of the Bg meson into an intermediate (virtual) B*. The 
Feynman diagram is shown in Fig.3. The subsequent process B* — » uu is determined in SM 
by the effective Hamiltonian eqn.(^, while the BqB*'y transition is taken as 

/: = ^?e.Aa/3e>X4*+h.c.. (10) 

This corresponds to the first part of eqn.(|^), which is for the CP-odd part of Bg. Hereafter, 
we will concentrate only on the CP-odd combination of Bg and Bg. The effective coupling 
constant is usually explained as a function of the invariant mass of the intermediate B*. 



Here we simply derive its value from the constituent quark model as a constant [11 



(11) 

3 \mb rugj 

The l/rrib term can also be neglected compared to the l/m^ term. Note that the potential 
model gives 

9--U^ + i-]^ (12) 



3 \Ab Aq, 



with 



(13) 



Arf = 0.59GeV, A;, = 4.93GeV (for Bj) 
A, = 0.66GeV, A^ = 4.98GeV (for Bs 
which is quite close to (|n]). Now the amplitude of the decay Bg •yw is: 

Here again we give the differential cross section for the purpose of being used experimentally: 



dV _ 2C^afl.ml, {m\ ^^^^ 



(im^p 3(487r)^m^^mg (^b* ~ ^^p)^ 
After integrating over the phase space, the decay width is obtained as 



(15) 
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where 



f{y) = -5y' + 62/' + I2y - 12 - 6(4 - y){l - y fHl - y). 



The formula (|T5i[T6|) are similar to the constituent quark model case eqn.(||,0). It is easy 



to see that if m^* = m^^, /b* = /s, are taken in eqn.(|T3|,|TB|), they will reduce exactly to 
eqn.(|,0). 

Now we use rriB* = 5.42 GeV and get: 

TiB, ^ ^uu) = 8.8 X 10-21 X GeV. (17) 

It is only slightly different from eqn.(D, because ttlb* is not quite different from rriB^- The 
branching ratio is 1.8 x 10^^, if /b* = fss = 0.2 GeV is used. The differential decay rate 
calculated in pole model is also displayed in Fig. 2 (dashed line) as function of m^p. 

For Bd meson, we get 

Y(B, ^ ^uu) = 1.1 X 10-21 X GeV. (18) 

The branching ratio is obtained as 2.4 x 10-^, also quite close to that obtained in constituent 
quark model. 

4 Conclusion 

We predict the branching ratios in SM for Bg —>■ '-fuu to be 10-^ and for B^ ^vv to be 
10~^. With these branching ratios, they are hopeful to be detected at future B factories or 
LHC. They can provide alternate channels for measuring Jb^ ifsj- 
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Figure Captions 



Fig.l Feynman diagrams in standard model for b — > quu. 

Fig. 2 Differential decay rate of Bg — >■ 'yi'i' versus the uu invariant mass m^p, with the 
solid line denoting results calculated from constituent quark model, and the dashed line 
corresponding to that of pole model. 

Fig. 3 Feynman diagram in pole model for Bq — > 'yi/i/. 
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